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ABSTRACT 

We present time-resolved optical photometry of the binary millisecond ‘redback’ pul¬ 
sar PSR J1023+0038 (=AY Sex) during its low-mass X-ray binary phase. The light curves 
taken between 2014 January and April show an underlying sinusoidal modulation due to the 
irradiated secondary star and accretion disc. We also observe superimposed rapid flaring on 
time-scales as short as ~20 s with amplitudes of ~ 0.1 — 0.5 mag and additional large flare 
events on time-scales of ~ 5 — 60 min with amplitudes of ^0.5-1.0 mag. The power density 
spectrum of the optical flare light curves is dominated by a red-noise component, typical of 
aperiodic activity in X-ray binaries. Simultaneous X-ray and UV observations by the Swift 
satellite reveal strong correlations that are consistent with X-ray reprocessing of the UV light, 
most likely in the outer regions of the accretion disc. On some nights we also observe sharp- 
edged, rectangular, flat-bottomed dips randomly distributed in orbital phase, with a median 
duration of ~250 s and a median ingress/egress time of ~ 20 s. These rectangular dips are 
similar to the mode-switching behaviour between disc ‘active’ and ‘passive’ luminosity states, 
observed in the X-ray light curves of other redback millisecond pulsars. This is the first time 
that the optical analogue of the X-ray mode-switching has been observed. The properties of 
the passive and active state light curves can be explained in terms of clumpy accretion from a 
trapped inner accretion disc near the corotation radius, resulting in rectangular, flat-bottomed 
optical and X-ray light curves. 

Key words: binaries: close - stars: fundamental parameters - stars: individual: 
PSR J1023+0038 - stars: neutron - X-rays: binaries 


1 INTRODUCTION 

The Galactic source PSR J1023+0038 (=AY Sex) was discovered 
by [Bond et all d2002t ). It was initially classified as a magnetic 
cataclysmic variable, but subsequent radio studies revealed the 
true nature of the compact object was revealed, with the dis¬ 
covery of 1.69 ms radio pulsations, establishing without a doubt 
that the comp act objec t was a rotation-powered millisecond pulsar 
dArchibald et alj2009t) . 

Millisecond pulsars are formed in X-ray binaries_ by ac¬ 
creting gas from a low-mass companion star ( lAlpar et all 
ll982|). The discovery o f accr et ing millisecond X-ray pulsars 
dWiinands & van der Klisl 1 19981 : IPatruno & Wattsl 120121) repre¬ 
sented a confirmation of this ‘recycling’ scenario for the long- 
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suspected evolutionary connection between low-mass X-ray bina¬ 
ries and rotation-powered recycled millisecond pulsars. The ac¬ 
creted matter spins up the neutron star to millisecond periods, but 
how and when the accretion process turns off remains unknown 
(see lTaurisl2012l ). PSR J1023+0038 has been heralded as the ‘miss¬ 
ing link’ between radio millisecond pulsars and low-mass X-ray 
binaries and has given new insights into this evolutionary process. 

In 2001 strong evidence for the presence of an accre¬ 
tion disc in PSR J1023+0038 was provided by optical observa¬ 
tions of double-peaked emission lines, accompanied b y short- 
time- s cale flickering and a blue optical spectrum dBond et al.1 
120021 : ISzkodv etak 120031) . From 2002 onwards optical/X-ray 
observations did not show the presence of an accretion disc 
dArchibald et alJl2009l . l20ld: iBogdanov etallkoi lh . Therefore, in 
2001, PSRJ1023+0038 was an accreting neutron star in a low- 
mass X-ray binary. Optical spectra taken in 2003/2004 showed 
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Figure 1 . Observed optical light curves of PSR J1023+0038 taken with the IAC80, INT and TNT (a-c) and the WHT (d). The light curves have been phase- 
folded using the ephemeris given in Section^ phase 0.0 is defined as inferior conjunction of the secondary star. The day/month and filter are also given in 
each panel. The solid line is a sinusoidal modulation (the secondary star’s ellipsoidal modulation plus heated star and accretion disc) scaled to fit the lower 
envelope (see Section^). © 2 Q14 RAS, MNRAS 000, ??-?? 
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Figure 2. The de-reddened and de-trended optical light curves of PSR J1023+0038, obtained by subtracting a fit to the lower-envelope of the light curves 
shown in Fig. |T| (a-c) and (d) show the data taken with the IAC80, INT, TNT and the WHT respectively. The solid circle marks the typical uncertainty in the 
data. The date day/month and filter are also given in each panel. 
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Table 1. Log of optical observations. The mean observed magnitude is given in the final column. The difference between the exposure time and time resolution 
is due the readout time of the detector. 


Telescope 

UT Date 

Start HJD-2450000 

Duration (h) 

Exp time (s) 

Time res. (s) 

Filter 

Seeing " 

Mag 

TNT 

2014/01/29 

6687.18760 

1.71 

0.97 

0.97 

r’ 

1.2 

16.63 

TNT 

2014/01/30 

6688.28886 

4.51 

0.97 

0.97 

r’ 

1.3 

16.16 

I AC 80 

2014/02/02 

6691.66953 

2.53 

15 

30 

r’ 

1.4 

16.31 

IAC80 

2014/02/03 

6692.55923 

5.02 

15 

30 

r’ 

1.5 

16.38 

IAC80 

2014/02/04 

6693.54567 

5.24 

15 

30 

r’ 

1.4 

16.61 

I AC 80 

2014/02/06 

6695.58626 

4.15 

20 

47 

r’ 

1.9 

16.79 

IAC80 

2014/03/01 

6718.39729 

7.43 

15 

30 

r’ 

1.7 

16.77 

IAC80 

2014/03/06 

6723.39929 

6.06 

30 

149 

r’ 

1.7 

16.80 

IAC80 

2014/03/08 

6725.39641 

5.37 

30 

45 

r’ 

1.9 

16.56 

WHT 

2014/03/16 

6733.35428 

0.88 

0.28 

0.31 

r’ 

1.3 

16.39 

WHT 

2014/03/16 

6733.35428 

0.88 

0.28 

0.31 

g’ 

1.5 

17.40 

WHT 

2014/03/16 

6733.35428 

0.88 

1.53 

1.55 

u’ 

1.7 

17.12 

IAC80 

2014/03/16 

6733.36879 

6.80 

30 

45 

R 

1.6 

16.21 

I AC 80 

2014/03/17 

6734.38586 

6.34 

30 

45 

R 

1.8 

16.66 

TNT 

2014/03/25 

6742.06245 

6.19 

3.37 

3.37 

r’ 

1.0 

16.60 

I AC 80 

2014/03/25 

6742.35805 

4.15 

30 

45 

R 

2.7 

16.52 

IAC80 

2014/03/26 

6743.34121 

6.60 

30 

45 

R 

1.8 

16.34 

TNT 

2014/03/28 

6745.10920 

4.98 

0.97 

0.97 

r’ 

1.1 

16.26 

I AC 80 

2014/03/29 

6746.35337 

6.37 

30 

45 

R 

1.6 

16.43 

INT 

2014/04/08 

6756.35347 

5.39 

10 

39 

r’ 

2.7 

16.61 

IAC80 

2014/04/24 

6772.37623 

4.08 

30 

45 

r’ 

1.3 

16.43 

WHT 

2014/04/29 

6777.35771 

2.17 

0.28 

0.31 

r’ 

1.0 

16.40 

WHT 

2014/04/29 

6777.35771 

2.17 

0.28 

0.31 

g’ 

1.0 

17.15 

WHT 

2014/04/29 

6777.35771 

2.17 

1.53 

1.55 

u’ 

1.1 

17.40 


only the absorption lines from the secondary star with no de¬ 
tectable emission lines and the radial velocity of the late-type G5 
companion star was found to be modula ted on the 4.754 h or¬ 
bital period f_P n ,.h: lThorstensen & ArmstrontJljQOsi) . Optical light 
curves taken in 2004 revealed a single humped modulation which 
can be explained in terms of an X-ray heated optical light curve 
dThorstensen & Armstrong 2005j). Between 2001 and 2007 radio 
pulsations were found ^Archibald et aljl2009l ) and until 2013 it has 
been consistently observed as an eclipsing radio millisecond pul¬ 
sar, where the radio eclipses are attributed to io nized material in 
the companion’s magnetosphere (Archibald et al. 2013). Given the 
~ 0.2 Mg se condary star, PSR J1023+0038 is known as a ‘red- 
back’ system jRobertsll201ll) . with a lower X-ray luminosity in the 
‘pulsar state’ due to a pu lsar-wind-driven shock near the inner La- 
grangian point feogdanov et alj201ll) . 

Since 2013 June, PSR J1023+0038 switched off as a radio pul¬ 
sar dStappers et al.ll20l3l : |Patruno et al.ll2014l) and when solar con¬ 
straints were overcome, optical spectroscopy in mid-October re¬ 
vealed strong double-pe aked Ha emission indicating that an accre - 
tion disc had re-formed dHalpem et al.ll2013l ; lLinares et alj|2014al) . 
Moreover, observations showed that the X-ray emission had in- 
creased by a factor of at least 20 compared to previous quiescent 
values |Kon^20D^_and the UV emission had brightened by 4 mag¬ 
nitudes (Patruno et al ■ 2 ml), with no detection of radio pulsations. 
All these factors provide compelling evidence for the presence of 
a recently formed accretion disc, indicating that PSR J1023+0038 
has switched back from a radio millisecond pulsar to a low-mass 
X-ray binary. The X-ray light curves during the accretion phase 
show rectangular flat-bottomed dips dTendulkar et al.l2014l) that are 
similar to the mode-switching behaviour between disc ‘active’ and 
‘pass ive’ luminosit y states observed i n othe r redback binary pul¬ 
sars dLinaresll2014l) . ICoti Zelati et al.l d2014l) determined the near- 
infrared to the X-ray spectral energy distribution (SED) and find 


that it is well represented by a model consisting of an irradiated 
companion, an accretion disc and intra-binary shock emission. The 
authors suggest a scenario in terms of an engulfed radio pulsar 
which is still active, but is undetectable at radio wavelengths due 
to a large amount of ionized material surrounding the compact ob¬ 
ject, material which is however not enough to totally quench the 
radio pulsar. 

In this paper we report on long-term and high time-resolution 
optical photometry obtained during the accretion disc/radio-quiet 
phase of PSR J1023+0038 in 2014. In a following paper we will 
prese nt the result of the spectroscopic campaign dLinares et al.1 

I2015I) . 


2 OPTICAL OBSERVATIONS AND DATA REDUCTION 

A total of 24 time-resolved optical light curves of PSR J1023+0038 
were obtained using four telescopes; the 0.8-m IAC80 telescope 
(Tenerife, Spain), the 2.4-m Thai National Telescope (TNT; Thai 
National Observatory, Thailand), the 2.5-m Isaac Newton Tele¬ 
scope (INT; La Palma, Spain) and the 4.2-m William Herschel Tele¬ 
scope (WHT; La Palma, Spain), during the period 2014 January to 
April. Conventional CCD optical light curves of PSR J1023+0038 
were obtained with the IAC80+CAMELOT and the INT+WFC 
with plate scales of 0.30 and 0.33 "/pixel, respectively. The ob¬ 
serving conditions were variable with seeing ranging from 1 to 
2". The exposure times used ranged from 10 to 60s with the 
r' or R filter depending on availability. For the WFC data the 
CCD was windowed to decrease the read-out time. Bias images 
and flat-fields were also taken. We used IRAI0 for our data re- 

1 IRAF is distributed by the National Optical Astronomy Observatory, 
which is operated by the Association of Universities for Research in Astron- 
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duction, which included b ias s ubtraction imag es and flat-fielding. 
The VAOPHOT software foeeg & Dovldl2001 1 was used to obtain 
differential light curves for PSR J1023+0038 and several compari¬ 
son stars by extracting the counts using variable aperture photome¬ 
try which scaled with the seeing. High time-resolution light curves 
were obtained at the TNT with ULTRASPEC j jDhillon_et al.l2014l) 
and at the WHT with ULTRACAM fohillon et al.ll2007lt . ULTRA- 


SPEC has an electron multiplier frame-transfer CCD with a pixel 
scale of 0.45 "/pixel and allows for single-band, high-speed imag¬ 
ing. ULTRACAM has three frame-transfer CCDs with a pixel scale 
of 0.3 "/pixel, allowing one to obtain simultaneous blue, green and 
red-band high-speed images. Due to the frame-transfer architecture 
of the CCDs, the dead-time is essentially zero with ULTRACAM 
and ULTRASPEC. With ULTRASPEC we obtained r'-band im¬ 
ages with exposure times of 0.96 and 3.35 s, with a typical photo¬ 
metric accuracy of 4.9 and 6.7 per cent, respectively. With ULTRA¬ 
CAM we obtained simultaneous it', g' and r'-band images with an 
exposure time of 0.28 s in g and r' and 1.53 s in u ', with a typi¬ 
cal photometric accuracy of 6.0, 7.0 and 5.1 per cent respectively. 
The ULTRACAM pipeline reduction software was used to debias 
and flat-field the data. The same pipeline was also used to obtain 
light curves for PSR J1023+0038 and several comparison stars by 
extracting the counts using aperture photometry with variable aper¬ 
tures that scaled with the seeing. 

For all of the data, calibrated magnitudes of the target with 
respect to a bright local standard SDSSJ102343.30+003819.1 
0=16.810, </=15.676, r'=14.773 and i'=14.952) were then de¬ 
termined. We checked that the local standard was not variable by 
comparing it to other bright stars in the field. For the data taken 
with the R filter, we used the SDSS transformations 0 to trans¬ 
form the r' magnitude of the local standard to R. We checked the 
transformation by using the IAC80 R and WHT r' light curves 
taken on March 26, which overlapped in time. We binned the light 
curves to the same time resolution and found both light curves to 
match with an offset of less than 2 per cent. As a check of the 
photometry and systematics in the reduction procedure, we also 
extracted light curves of a comparison star of similar brightness 
to PSR J1023+0038, which showed no significant variability above 
what is expected from photon statistics. 


3 OPTICAL LIGHT CURVES 


We determined the colour excess by using the empirical relation 
between the colour excess and the t otal column density of hy¬ 
drogen Nh IPredehl & Schmitt] HWSjj-Using the_observed value 
for Nh = 5.2 x 10 20 cm ~ 2 ICoti Zelati et all 1 20141) a nd the ra¬ 
tio Av/E(B — U)=3.1 i lCardelli. Clavton. & Mathisl 1 19891) we 
obtain a colour excess of E(B-V)=0.073. Fig.Q] shows the de- 
reddened flux IAC80, INT, TNT and WHT light curves where 
we have used the colour excess determined above and then con¬ 
verted the magnitudes to flux density using the appropriate ze- 
ropoint for each band. We have also phase folded the data us¬ 
ing t he ephemeris P or b=0. 1980963569 d and To =2456642.6335 
HJDl Archibald et aPl2009l : iLinares et alJl20l3) . where phase 0.0 
is defined as inferior conjunction of the secondary star. 


omy, Inc., under cooperative agreement with the National Science Founda¬ 
tion. http://iraf.noao.edu/ 

2 http://www.sdss.org/dr5/algorithms/sdssUBVRITransform.html 
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Figure 3. The February 6 (bottom), March 25 (middle) and April 8 (top) 
r'-band de-reddened and de-trended light curves. The March 25 light 
curve taken with TNT+ULTRASPEC has been re-binned to the same time- 
resolution as the IAC80 and INT light curves. The variability represents 
genuine variations of the target as indicated by the error bar plotted in the 
top right hand corner of each panel. 


3.1 Flaring activity 

The light curves show dramatic variations from night to night su¬ 
perimposed on an underlying sinusoidal modulation, which is due 
to a combination of the second ary star's e llipsoidal modulation and 
it’s heated inner face jThorstensen & Armstrondl2005h . On some 
nights there is additional rapid flaring (February 6; March 25; April 
8) and on other nights there are further additional large flare events 
(January 1, February 4, March 8). In order to isolate the flaring ac¬ 
tivity (by removing the shape of the secondary star's modulation) 
we subtract a scaled version of the underlying sinusoidal modula¬ 
tion, defined by fitting the lower-envelope of the light curve taken 
on March 17 with a sine-wave. The March 17 light curve was cho¬ 
sen because it clearly showed the sinusoidal modulation with lit¬ 
tle flaring behaviour. Fig.[2] shows the de-reddened and de-trended 
IAC80, INT, TNT and WHT flare light curves. 

Fig.[3] shows the flares on February 6, March 25 and April 8. 
The variability represents genuine variations of the target. Notice¬ 
able from the IAC80 and INT data are numerous rapid flare events 
that are not resolved. The TNT March 25 data were taken with 
ULTRASPEC and so have a much higher time-resolution than the 
IAC80 and INT data (see TableQ} and so the flares are better re¬ 
solved. In order to directly compare the light curves, the March 25 
data were re-binned to the same time-resolution as the other data. 
The subsequent light curves obtained are very similar and hence 
we can conclude that the IAC80 and INT data do not resolve the 
rapid flares. The TNT data show rapid flare events on time-scales of 
^minutes with amplitudes of ~ 0.1 — 0.5 mag (see Fig.[4]l. Large 
flare events on time-scales of ~ 5 — 60 min with amplitudes of 
~ 0.5 — 1.0 mag are also seen (e.g. February 2, March 1) super¬ 
imposed on the rapid flare events. In Table[3]we give the fractional 
root-mean-square (RMS) of the variability, determined after bin¬ 
ning the light curves to the same time-resolution of 60 s. 
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Figure 4. (a-c) The de-reddened r'-band light curves of PSRJ1023+0038, taken with the TNT+ULTRASPEC and the WHT+ULTRACAM. The right panel 
in each plot shows a histogram of the flux values after correcting for the secondary star’s trend. The mode-switching active and passive state behaviour is only 
observed in the January 29, 30 and April 29 light curves. The vertical dashed lines shown marks the cut-off flux value that defines the active flare and active 
passive respectively, (d) Histograms of the observed rectangular mid-dip phase, ingress/egress times, duration and the separation between consecutive dips. 



Figure 5. Left: A close-up of the flat-bottomed dip features observed in the April 29 WHT light curve. We show the de-reddened r'and the flux ratio f u / / f g t 
and fgt / f r / light curves. Notice the increase in the blue colour during the dips. The error bar marks the typical uncertainty in the data. Right: Colour-magnitude 
diagram for the April 29 WHT data. We show the de-reddened r 'and the flux ratio f u / / f g / and f g / /f r / data. The solid black and red points show the active 
and passive state data respectively. Note the increase in colour during the dips, i.e. they become more blue. 
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3.2 Rectangular dips 

Fig-0shows the observed r'-band light curves and flux histograms 
taken at the highest time-resolution; 0.31 s for the WHT data and 
0.97 or 3.37 s for the TNT data. These light curves, taken on Jan¬ 
uary 29, 30 and April 29, show rectangular, flat-bottomed dip fea¬ 
tures which are similar to the mode-switching behaviour (passive 
state and active state) that have been observed in the X-ray light 
curves of PSR J1023+0038 dTendulkar et all2014l) and other X-ray 
binary pulsar (see lLinaresll2014l . and references therein). The his¬ 
togram of the flux values (using the flare light curves described in 
Section[3]with the mean de-reddened flux added to preserve the flux 
values) clearly show a bimodal distribution between the passive and 
active -states (see Fig.[4]). 

In total we identify 16 rectangular dips in the January 29, 30 
and April 29 light curves. In order to estimate the properties of 
these dips, we fit each feature with a rectangular function, allow¬ 
ing us to determine the ingress, egress, duration and centre of the 
dip segment. We find that the dips are generally symmetrical with 
ingress and egress times in the range 12 to 35 s (median of ~ 20 s), 
a dip duration in the range 80 to 1300 s (median of ~ 250 s) and a 
dip separation in the range 200 to 1900 s (median of ~ 900 s). We 
also compute histograms of the observed dip phase, the duration of 
the dips and the separation between consecutive dips (see Fig|4j, 
however, the low number of dips observed implies that we cannot 
comment on the statistical properties of the dips with confidence. 
In Fig.[5]we show a zoom of a typical rectangular dip feature ob¬ 
served on April 29 with ULTRACAM. The feature is flat-bottomed 
with an increase in the /„/ /f g i and f g i /f r i colour ratios during 
the dip, i.e. the dips become more blue. We also show the colour- 
magnitude diagram of the whole light curve. The increase in colour 
ratio during the dips is clearly seen. 


4 X-RAY/ULTRAVIOLET LIGHT CURVES 


We analyzed all Swift observations of PSR J1023+0038 in the on¬ 
going accretion phase, taken between 2013 October 18 and 2014 
June 11, including data from both the X-ray (XR T, iBurrows et al.l 
120051) and UV/optical (UVOT. lRoming et al.l2005h telescopes. This 
amounts to 52 observations, an accumulated XRT exposure of 
94 ks and a total of 84 UVOT images acquired in four differ¬ 
ent filters. Smaller sets of Swift observations of PSR J1023+0038 
have been analyzed elsewhere (|PatrunoetalJ 120141 : iTakata et al.1 
l2014l: iLinarejEobl ICoti Zelati et al .1120 1 4) . We used the quasi- 
simultaneous XRT/UVOT average flux per observation (which lasts 
typically between 1 and 3 ks), as we are interested in the long-term 
(weeks-months) X-ray-UV correlation. Because the UVOT filter 
configuration and total XRT exposure vary, each observation had 
UVOT data collected in 1 to 3 different filters. 

All XRT observations were taken in photon counting (PC) 
mode. For the three brightest observations (on 2014 February 27, 
March 3 and 11, during which the count rate reaches ~0.6 cts s -1 ), 
we used a 20-pixel radius extraction region and excluded the inner¬ 
most pixels (within a 2-pixel radius) to avoid pile-up. The count 
rate during the remaining observations is <0.5 cts s -1 (below 
the PC pile-up threshold) and we used a 30 arcsec radius circular 
extraction region. We subtracted the background using a nearby 
60 arcsec radius source-free region, and grouped the source spec¬ 
tra to a minimum of 15 counts per channel. We created ancilliary 
response files for each observation, applying the vignetting correc¬ 
tion and using the corresponding exposure maps. After incorporat¬ 
ing the latest XRT response matrices, we fitted each spectrum in the 


Table 2. Properties of the simultaneous X-ray/UV observations of 
PSR J1023+0038. N is the number of simultaneous X-ray/UV observa¬ 
tions, R is the correlation coefficient and /3 is the power-law index of the 
form Lv } uv oc L^. 


Band 

N 

R 

P 

U 

1 

0.26 

0.04(4) 

UVW 1 

35 

0.62 

0.77(2) 

uvm2 

14 

0.89 

1.13(4) 

uvw2 

22 

0.87 

1.02(3) 


0.5-10 keV band using XSPEC (v. 12.7.1 . lArnaucfaoOfil . except two 
spectra which had less than 45 source counts and were not fitted). 
We used a pow e r-law modelcorrected for absorption (with TBABS, 
Iwilms. Allen. & M cCray! 200Ch . fixing the column density at the 
value measured bv ICoti Zelati et aD J2014 ) from the average XRT 
spectrum (7Vh=5.2x 10 20 cm -2 ). We calculated the 0.5-10 keV X- 
ray luminosity from the unabsorbed flux in the same band using the 
distance of 1.37 kpc JPeller et alj2012i) . 

After checking that the aspect correction was applied, we 
summed all extensions within each UVOT image mode file (using 
the Swift ftool UVOTIMSUM, v. 1.6) to obtain one image per filter 
and observation. Because our ULTRACAM and ULTRASPEC data 
provide much better measurements of the fast optical variability of 
PSR J1023+0038 (see Section[3}, we do not include in thisjmaly- 
sis UVOT high-time resolution (event mode) data (ICoti Zelati et al.1 
l2014h . We then performed aperture photometry on each summed 
image using the standard 5 arcsec aperture radius, and a nearby 
source-free background region of 12 arcsec radius. 

In Fig.[6]we show the Swift XRT 0.5-10 keV unabsorbed long¬ 
term X-ray light curve. The X-ray light curve shows considerable 
activity with an average luminosity of Lx ~ 2 x 10 33 erg s -1 . To 
search for possible X-ray/UV correlations in Fig.[6]we determine 
the simultaneous unabsorbed X-ray 0.5-10 keV and de-reddened 
U, uvw 1, uvm2 and uvw2-band observations. The UV data were 
de-reddened using the colour excess derived in Section[3] The X- 
ray and uvw 1, uvm2 and uvvtZ-band data are clearly correlated, with 
a significance of >99 per cent. However, given that there are only 7 
points in the X-ray- U correlation, the correlation is only significant 
at the 50 per cent level. We fit the X-ray/UV data with a power-law 
of the form L„ ; uv oc L^ and determine the power-law index (see 
Table[2}. 


5 OPTICAL POWER DENSITY SPECTRUM 

We compute the power density spectrum (PDS) of the de-reddened 
and de-trended light curves (see Section[3] with the mean de- 
reddened flux added to preserve the flux values) using the Lomb- 
Scargle method to compute the periodograms ( iPress et al.l 119921) 
and the same normalization method as is commonly used in X- 
ray astronomy, where the power is normalized to the fractional 
RMS amplitude per Hertz dvan der Klislll989h . To compute the 
PDS we use the constraints imposed by the Nyquist frequency, 
th e typical duration of each observation and use the recipe given 
in iHorne & Baliunasl (11986h to calculate the number of indepen¬ 
dent freque ncies. We then bin and fit the PDS in logarithmic 
space JPapadakis & Lawrencelfl993 ) where the errors in each bin 
are determined from the standard deviation of the points within 
each bin. The white noise level is subtracted by fitting the high- 
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Figure 6. Swift XRT andUVOT light curves of PSR J1023+0038 between 2013 Oct 18 and 2014 June 11. Left: Swift XRT 0.5-10 keV light curve. The arrows 
mark the time of our IAC80, INT, TNT and WHT optical data. Right: Swift simultaneous UVOT (U, uvwl, uvw2 and mm2) and XRT 0.5-10 keV observations. 
The error bar plotted in the top left hand comer of each panel. Tthe uncertainty in the UVOT data is smaller than the symbol size. A power-law fit to the data 
of the form uv oc /,{! is shown. 


Table 3. The fractional RMS of the r'-band flare light curves and power- 
law (PL) fits to the power density spectra. The light curves where active and 
passive state are seen are not included. The mean de-reddened flux values 
are given and the fractional RMS values have been determined after binning 
the light curves to the same time-resolution of 60 s. 


Telescope 

UT Date 

Flux (mJy) 

RMS (%) 

PL index 

TNT 

2014/01/29 

0.94 

9.5 

-1.43± 0.04 

TNT 

2014/01/30 

1.53 

30.0 

-1.54± 0.04 

IAC80 

2014/02/02 

1.29 

30.1 

-1.53± 0.07 

IAC80 

2014/02/03 

1.13 

21.7 

-1.40± 0.04 

IAC80 

2014/02/04 

0.91 

21.1 

-1.39± 0.04 

IAC80 

2014/02/06 

0.82 

10.3 

-0.61± 0.07 

IAC80 

2014/03/01 

0.79 

21.0 

-0.97± 0.03 

IAC80 

2014/03/06 

0.80 

6.6 

-0.91 ±0.08 

IAC80 

2014/03/08 

1.00 

27.1 

-1.34± 0.05 

IAC80 

2014/03/16 

1.23 

34.6 

-1.29± 0.05 

IAC80 

2014/03/17 

0.81 

8.0 

-0.65± 0.05 

TNT 

2014/03/25 

1.53 

29.9 

-1.05± 0.01 

IAC80 

2014/03/26 

0.86 

8.2 

-1.22± 0.05 

TNT 

2014/03/28 

1.32 

23.8 

-1.83± 0.01 

IAC80 

2014/03/29 

1.35 

15.0 

-1.14± 0.05 

INT 

2014/04/08 

1.24 

9.3 

-0.72± 0.05 

IAC80 

2014/04/24 

1.11 

9.9 

-1.06± 0.07 


est frequencies (>300 Hz) with a white noise plus red noise model. 
As an example of the typical PDS, in Fig.[7Ji we show the PDS 
of the light curve taken on March 8 and in Table[3]we give the 
power-law index obtained for each light curve. We find that the 
PDS of the flare light curves is dominated by a red-noise com¬ 
ponent with a median power-law index of -1.2, typical of ape¬ 
riodic activity in X-ray binaries and X-ray transients in outburst 
and quiescence, which have a power-law index in the range - 


1.0 to -2.0 ( 

Zurita, Casares, & Shahbaz 2003; Shahbaz et al. 2003 

Hynes et al 

2003ajb; Shahbaz et al. 20041. 2005, 20ll 201J. 


In order to determine the significance of possible peaks above 
the red-noise level^we use a Monte Carlo simulation similar to 
IShahbaz et al.1 (2005 ). We generate light curves with exactly the 
same sampling and integration times as the real data. We start 
with a model noi s e ligh t curve generated using the method of 
iTimmer & Koenigj 119951 ). with a power-law index as determined 
from the PDS of the observed data and then add Gaussian noise us¬ 
ing the errors derived from the photometry. We compute 1000 sim¬ 
ulated light curves and then calculated the 99 per cent confidence 
level at each freque ncy taking int o account a realistic number of 
independent trials l lVaughanl 120051) . As an example, in Fig.[7ji we 
show the 99 per cent confidence level contour, which rules out any 
significant peaks in the PDS. 

We also compute the PDS of the multi-band ULTRACAM 
data taken on April 29, which shows strong active/passive state be¬ 
haviour. We find the power-law index of the PDS to be -1.43±0.05, 
-1.34±0.07 and -1.10±0.05 in the r', g 'and it'bands respectively, 
which implies that there is more high frequency variability at 
shorter wavelengths. 

For the light curves that show mode-switching behaviour (Jan¬ 
uary 29, 30 and April 29), we compute the PDS of light curves 
during the active and passive states to see if they are different. We 
use the flare light curves computed in Section[3]with the mean de- 
reddened flux added to preserve the flux values. In order to iso¬ 
late the passsive- and active states, we determine the cut-off using 
the histogram of the r'-band flux values, where the cut-off is the 
flux value that divides the bimodal distribution (see Fig.[4]>. For the 
WHT data the times in the r'-band light curve corresponding to 
flux values greater than or less than the cut-off value allows us to 
define the time values in the g and u'-band active and passive state 
light curves. Although we find a number of interesting differences 
in the slope of the active and passive state PDS, a Monte Carlo 
analysis shows that the differences are due to the window function 
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(a) IAC80 March 8 


(b) 


TNT January 29 


TNT January 30 


o 

a_ 




(C) 


WHT April 29 


WHT April 29 





Figure 7. (a) An example of the typical power density spectrum (PDS) of the IAC80 and INT light curves. The PDS shown is for the r'-band light curve 
taken on March 8 and the solid line represents the 99 per cent confidence level, which rules out any significant peaks in the PDS. (b) The PDS of the r'-band 
TNT light curves taken on January 29 and 30. (c) The PDS of the r' and g' and u'-band WHT light curves taken on April 29. In each plot the dot-dashed lines 
show a power-law fit to the PDS. 


of the light curves. We perform a simulation in which we simu¬ 
lated a red-noise light curve with the same time values, power-law 
index, mean and rms as the r' band light curve. We then selected 
active and passive sections (the same as what we observe) and the 
compute the power spectrum and determine the power-law index. 
These were then compared to the power-law index of the simulated 
r'band PDS. One would expect the PDS of the active and passive 
light curves to be the same as the PDS of the whole light curve. 
However, we find that the power-law index of the active and pas¬ 
sive light curves are very similar, but different to the power-law 
index of the whole light curve. 


6 PROPERTIES OF THE LARGE FLARE 

The multi-colour data taken with ULTRACAM allow us to deter¬ 
mine the colour of the large flare event and dips (see Fig.[2ji). The 
March 16 light curve starts off with relative low variability (ac¬ 
tive state) which develops into a large flare event lasting at least ~ 
30 min. In contrast the April 29 light curve shows dips in flux. To 
determine the colour of the flare event we use the flare light curves 
computed in Section[3] where the secondary star’s sinusoidal mod¬ 
ulation is removed and the mean de-reddened flux value is added to 


preserve the flux values in each band. For the March 16 data we cal¬ 
culate the average active state and near-peak flux values for the r', 
g 'and «'-band light curves. The flux of the flare is then obtained by 
subtracting the active state flux from the near-peak flux value. For 
the April 29 data we calculate the average active and passive state 
flux values. The flare of the dips is then obtained by subtracting the 
active state flux from the passive state flux (see Table[4}- 

In an attempt to interpret the broad-band spectral properties 
of the flare event, we compare the de-reddened colours with the 
prediction for different emission mechanisms, such as a blackbody, 
an optically-thin layer of hydrogen and synchrotron emission. The 
most likely model for a thermal flare is emission from an optically 
thin layer of recombining hydrogen, which is the mechanism gener¬ 
ally accep ted for stellar flares in single stars and Cataclysmic Vari¬ 
ables (see lKerr & Fletcherll20 1 4l and references within). We there¬ 
fore determine the continuum emission spectrum of an LTE slab of 
hydrogen for different baryon densities, (Vh and temperatures Th. 
For each model we then compute the given emission spectrum and 
the expected flux density ratios f u i / f g i and f g > /f r r using the syn¬ 
thetic photometry package SYNPHOT (IRAF/STSDAS). Given the 
intrinsic model flux we can then determine the corresponding ra¬ 
dius ( Ra ) of the region that produces the observed de-reddened 
flux at a given distance of 1.368 kpc foeller et alj|2012l ). 
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25 24 23 22 log 10 N H 



Figure 8. The colour-colour diagram of the WHT light curves of 
PSRJ1023+0038 taken on the March 16 (blue) and April 29 (green). The 
colour during the active state is represented by a solid square, whereas the 
colour of the large flare event and dips is represented by a solid circle. 
The dashed lines show optically-thin hydrogen slab models for different 
baryon column densities log ,Yh , where the numbers mark the temperature 
in 100 K units. From right to left we show models with log Ah =22 — 
30 nucleons cm -2 . The solid line shows a power-law model ( F „ oc i/*) 
with indices a=—2.0 to —4.0 marked as stars. The dotted line shows a black- 
body model marked in 100 K units as open squares. The large red star shows 
the colour of a G5 main sequence star. 


In Fig .[8] we show the colours of the March 16 active and 
flare-state data and the expected colours for different emission 
models. We can see that the blackbody or synchrotron emission 
models do not reproduce the observed flux ratio’s and so these 
models can be ruled out. For the March 16 light curve, the ac¬ 
tive state flux ratios (see TableO correspond to a region with 
Ah ~ 10 25 nucleons cm -2 and Th ~ 6600 K. The flux ratio 
of the flare (near-peak flux minus active state flux) corresponds 
to a region with Ah ~ 10 23 nucleons cm” 2 , Th ~ 8400 K and 
Rh ~ 0.35 Rq. Hence we find that during the large flare event 
there is a decrease in baryon density and an increase in tempera¬ 
ture, suggesting that the flare emission becomes ho tter and more 
optically thin. Using the distance of 1.368 kpc foeller et alj|2012h 
and the corresponding filter bandwidth, the optical peak luminosity 
of the large flare is 4.9, 2.9 and 3.3x 10 32 erg s” 1 in the u', g' and 
r' -bands respectively. Although we do not observe the full decay 
of the flare event, we estimate that the flare lasts ~ 1000 s (~ 0.06 
in phase) with an energy of 6x 10 35 ergs in the combined u , g' and 
r'-bands. 

For the April 29 light curve, the active state flux ratios cor¬ 
respond to a region with Ah ~ 10 25 nucleons cm -2 and Th ~ 
6900 K (similar to the active state in the March 16 light curve). 
The flux ratio of the dips (i,e. the flux ratio between bands after 
subtracting the active state flux from the passive state flux in each 
band) corresponds to a region with Ah ~ 10 3 ° nucleons cm” 2 , 
Th ~ 4600 K and Rh ~ 0.35 Rq. Note that if we assume no 
reddening the colour ratios change by at most 6 per cent. 


7 DISCUSSION 
7.1 The flaring activity 

IPatruno et alj l l2014l) reported the X-ray and UV/optical signatures 
of a state change in PSRJ1023+0038 during mid 2013, which 
coincided with the disappearance of radio pulsations and the X- 
rays show a softer spectrum than during the millisecond pulsar 
state, indicating a change in the origin of the X-ray emission. 17- 
band variability of ~ 0.3 mag on time-scales of ~ 1 h was seen as 
well as X-ray flux changes by up to a factor of ~100 on time- 
scales less than < 10 s. During the previous active phase in 2002, 
slow (hours) and rapid^minutes to seconds) optical flickering was 
observed feond et alj 12002 ). Our light curves observed in 2014 
during the accretion-powered phase show rapid flaring events on 
time-scales of minutes with amplitudes of ~0.1-0.5 mag and/or 
large flare events on time-scales of ~5-60 min with amplitudes of 
~0.5-1.0mag (see Section [3Tt . The PDS of the flare light curves 
has a median power-law index of —1.2 (see Section[5), t ypical 
of X-ray binaries llZurita. Casares. &~Shahbaz 2003; Shahbaz et al] 

I 2003 I I 2004 L I 2005 LI 20 1 d. | 20 1 3l) . 

Since the accretion disc is sandwiched by a hot, optically thin, 
fully ionized corona, two possible origins for the flares are (1) the 
optically thick, geometrically thin accretion disc or (2) from the 
corona, leading to enhanced energy release over a short period of 
time. In the context of the propeller effect, the inner edge of the 
accretion disc rotates slower than the neutron star magnetosphere 
and in-falling matter is no longer able to accrete onto the neu- 
tron star, but instead is blo wn out of the system by a propeller 
(Illarionov & Su nvaevlll975h . leading to flaring activity. This has 
been used to explain the large luminosity varia tions in the outburst 
tails of SAX J1808.4-3658 fcampana. Stella. & Kenne~all2008l) . In¬ 
deed the sporadic large flaring events observed in PSR J1023+0038 
during its accretion-powered phase could also be explained in 
terms of a propeller effect, leading to the ejection of hot opti¬ 
cally thin fire balls, as is o bserved in some cataclysmic variables 
dPearson, Home. & Skidmorell2003l) . 

For the large flare event observed on March 16 (lasting at 
least ~ 30 min) where we have multi-colour information, we find 
that the event arises from a relatively optically thin region with ~ 
8700 K, Ah ~ 6x 10 22 nucleons cm” 2 and a radius of ~ 0.35 Rq; 
the out-of-flare state is cooler and more optically thick (see Fig.|8j. 
The energy of the flare in the optical bandpass of 6xl0 35 ergs 
(see Section^ is two orders of magnitude larger than the op¬ 
tical lumin osity of flares produce d by chromospheric activity in 
flare s tars dWalkowicz et al.ll201lh . Ichurch & B aiucinska-Churcl] 
(2004) have measured the radial extent of the accretion disc corona 
(ADC) in low mass X-ray binaries. They find the ADC extends 
from 20,000 km in the faintest sources to 700,000 km in the bright¬ 
est, typically 15 per cent of the accretion disc radius. In this context, 
our data suggest that the large flares arise from a large fraction of 
the accretion disc, most likely from the accretion disc corona. 


7.2 The X-ray/UV correlations 

Power-law correlations between the optical and X-ray luminosities 
in an X-ray binary are naturally expected. The optical emission in 
low-mass X-ray binaries is generally thought to arise in the outer 
accretion disc as the result of X-ray reprocessing and the luminosi¬ 
ties are correl ated in the f orm T „ : y oc (V refers to the V- 
band) dvan Paradiis & McClintockll 1 994l ). The index (3 depends on 
the surface brightness of the accretion which is given by S v oc T“, 


© 2014 RAS, MNRAS 000, ??-?? 






























Optical observations ofPSR J1023+0038 11 


Table 4. Mean flux values of the WHT light curves taken on March 1 6 and April 29. For the March 1 6 event the mean active, near peak fluxes are given. The 
flare flux is calculated by subtracting the near-peak flare flux from the mean active state flux and we also give the optical luminosity of the flare. For the April 
29 light curves the mean active and passive state fluxes are given as well as the X-ray (0.5—lOkeV) to op tica l (u' , g' and r') luminosity ratio, where the X-ray 
active and passive state luminosities are 2.85xl0 33 ergs -1 and 5.70xl0 32 ergs -1 respectively dl.inares|[2014l) 


UT Date 

Band 

Active state 

Near-peak 

Peak—Active 

Z/OPT 



flux (mJy) 

flux (mJy) 

flux (mJy) 

(ergs -1 ) 

March 16 

fu' 

0.515(2) 

1.369(7 

0.855(9) 

4.9xl0 32 


h' 

0.398(1) 

0.868(4) 

0.471(4) 

2.9xl0 32 


Sr' 

0.973(2) 

1.921(8) 

0.948(8) 

3.3xl0 32 


fu'/fg' 

1.294(7) 

1.576(12) 

1.825(15) 



fg’/fr' 

0.409(1) 

0.452(3) 

0.497(2) 




Active state 

Passive state 

Active—Passive 

Lx/Lopt 



flux (mJy) 

flux (mJy) 

flux (mJy) 

Active Passive 

April 29 

fu' 

0.563(2) 

0.538(1) 

0.025(2) 

9.0 1.9 


fg' 

0.720(1) 

0.632(1) 

0.087(2) 

6.5 1.5 


Sr > 

1.375(2) 

1.174(1) 

0.202(3) 

6.2 1.4 


fu'/fg’ 

0.523(1) 

0.539(1) 

0.283(21) 



fg’ Hr' 

0.782(2) 

0.851(1) 

0.432(8) 



where P = a /4 Ivan Paradiis & McClintockll 1994 1. In the V-band 
a ~ 2.0 and so P ~ 0.5 Ivan Paradiis & McClintockll 19941 ). It can 
be shown that a increases with decreasing wavelengths; in the UV- 
band, a ~ 3.7 (/3 ~ 0.9), in the U-band a ~ 2. 7 (P ~ 0.7) and 
in the Tv-band a. ~ 1.2 (/3 ~ 0.3). lRussell et al.l 120061) show that 
similar correlations are expected when the optical flux originates in 
the viscously heated disc because both the X-ray and optical fluxes 
are linked through the mass accretion rate. For a viscously-heated 
steady-state disc, with outer disc temperatures of 8000-12,000 K, 
P is dependent on wavelength and increases at shorter wavelength; 
P ~ 0.3 in the 7\-band, p ~ 0.5 in 17-band and P ~ 0.60 
in the 17-band. Note that the power-law index due to an irradi¬ 
ated accretion disc is steeper than that due to a viscously heated 
disc, for a given waveband. iRussell et alj 120061) have quantified 
the disc and jet contribution from the optical/infrared and X-ray 
(2-10 keV) correlation in X-ray binaries in the hard state. From 
an observational point of view a strong optical/X-ray luminos¬ 
ity correlation from quasi-simultaneous data has been observed in 
neutron-star X-ray binaries in the hard state, which takes the form 
L^oir oc T,^ 63 * 0 - 04 and holds over seven orders of magnitude in 
Lx (where OIR covers the B to 77-band). The authors concluded 
that X-ray reprocessing dominates the optical/near-infrared emis¬ 
sion in neutron-star X-ray binaries in the hard state, with possible 
contributions from a jet (only at high luminosity) and a viscously 
heated disc. 

From Fig.[6] we find the UV and X-ray luminosities are 
strongly correlated with a power-law index ~ 1.0 in the UV bands 
(see Table0. This is steeper than what is expected from a viscously 
heated disc (0.6) and is more in line with the index expected from 
an irradiated accretion disc (0.9). 


7.3 The rectangular dips 

X-ray mode-switching behaviour observed as fast transitions be¬ 
tween passive and active s tates (rectangular dip s ), has been de- 
tected in PSRJ1023+003 8 ITendulkar et al. 20141 : [Archibald et al.l 
120141) and other similar systems iLinaresll2014l) . These dips, which 
are rectangular, are distinct from the dipping activity observed 
in some high inclination low-mass X-ray binaries, called X-ray 


dippers, where intense dipping activity at specific orbital phases 
are observed due to azimuthal absorbing material located close to 
the accretion disc bulge Iwhite, Nagase, & Parman 119951) . XMM- 
Newton observations of PSRJ 1824-24521, with an X-ray lumi¬ 
nosity 10 35 ergs -1 , showed rectangular dips with ingress and 
egress times of ~ 200 s and dip durations of up to a few thou¬ 
sand seconds iFerrigno et alj|2014l) . Chandra observations in qui¬ 
escence showed mode-switching behaviour with states las ting for 
~ 10 h and transition times of less than 500 s iLinares et al . 2014bh . 
Dips which show a spectral hardening, ingress and egress times 
of <10s and durations between 200 and 800 s have also been 
observed in XMM-Newton observations of XSS J1227.0-4859 
Ide Martino et al.ll2013l) . iLinaresI 120141) has shown that all three 
redbacks observed so far which have transitioned to the accre¬ 
tion phase (PSRJ1023+0038, XSSJ 1227.0-4859 and PSRJ182- 
24521) show X-ray mode-switching, with little or no change in 
X-ray photon index. For XSS J1227.0-4859 and M28-I, where 
there are measurements of Nu during the a ctive and passive states , 
no significant chang e in Vh is observed (de Martino et al.l [20131 : 
ILinares et alJl2014a) . NuSTAR observations of PSR J1023+0038 in 
October 2013, with an average X-ray luminosity of 5.8 x 10 33 
ergs -1 , show sharp-edged, rectangular, flat-bottomed dips with 
durations of 30 to 1000 s, ingress and egress times of 30 to 60 s 
and typical separations of ~ 400 s ITendulkar et aljfiouil) . The X- 
ray luminosity of the activ e and passive states only changes by a 
factor of ~5 lLinaresl2014l) . 

Our high time-resolution optical observations show sharp- 
edged, rectangular dips (see Section [3~2l i. which is the first time 
that the optical analogue of the X-ray active/passive state mode¬ 
switching has been observed. The histogram of the optical mid-dip 
position as a function of orbital phase is relatively flat (see Fig[4]l, 
suggesting that the dips are uniformly and randomly distributed 
with orbital phase, similar to what has been observed in the X-rays. 
The timescale properties (ingress/egress, separation etc) of the op¬ 
tical and X-ray dips are very similar, which suggests that we are 
seeing echoes of the X-ray variability. We observe an increase in 
the f u i / f g i and f u i / f g i colour ratio during the dips, which is due 
to the disappearance of a red spectral component (see FigO. The 
X-ray photon index of the active and passive state X-ray data does 
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not change significantly lLinares|[20l4) . The decrease in the X-ray 
luminosity between the active and passive states is a factor of 
dLinare sl2014h . much larger than the decrease in optical luminosity 
between the active and passive states; ~10 per cent; the decrease in 
the X-ray-to-optical luminosity ratio during the active and passive 
state of a facto r between 4.4 and 4.7 (see Tabled. 

[Bogdanov et all 120141) recently presented the results of a 
multi-wavelength study of PSR J1023+0038 during it’s low-mass 
X-ray binary state. B-band light curves taken with a cadence of 13 s 
during late 2013 and 2014 clearly show the secondary star’s heated 
ellipsoidal modulation with superimposed flaring activity, similar 
to our optical light curves (see Section[3j. The authors also obtained 
considerable simultaneous XMM-Newton X-ray and B-band light 
curves taken with a time resolution of 10 s. They find that none 
of the B-band light curves show the mode switching behaviour as 
seen in the X-rays. At first light this seems to contradict our find¬ 
ings, where we clearly observe optical mode switching in some of 
our high-time resolution light curves. However, it should be noted 
that in the ULTRACAM data, where mode switching is observed, 
the mode switching is strongest in the r'band and is significantly 
weaker in the g' and w'-bands; the mode switching flux changes by 
a factor of 17, 14 and 5 per cent in the r', g' and it'-band light curves 
respectively (see Fig.jTJi and Table[4j. Therefore, it is not surpris ing 
that the B-band light curves obtained bv iBogdanov et al] |2014|) do 
not show cle ar discernible mode switching behaviour. 

[Bogdanov et alj d2014l) extensive XMM-Newton observations 
allowed them to and observe the source in various states: ~20, ~70 
and ~10 percent of the time in the low- (passive), high- (active) 
and flare-state respectively. Our ground-based optical observations 
total about 96 hr. During this time, we find the source to be in the 
flare-state for ~ 10 hr. Although we identify passive/active state 
behaviour in only 6 hr of data, it should be noted that the poor time- 
resolution of the IAC80 and INT data does not allow us to resolve 
the rapid flaring. It is clear that more extensive high-time resolution 
optical data is needed before we can compare the properties of the 
optical and X-ray state changes. 


7.3.1 Proposed scenario for the optical and X-ray 
mode-switching 

The details of the interaction between the rotating magnetosphere 
of the ne utron star and the inner ac cretion flow are complicated 
(e.g., see IUzdenskvll2004l ; lLaill20l4 However, it is generally as¬ 
sumed that accretion onto the neutron star’s surface occurs when 
a geometrically thin disc ( truncated at the ‘magnetospheric ra¬ 
dius’ r m ; where the ram pres sure of the gas balances the pres¬ 
sure of the magnetic field, [Pringle & Reesll 19721 ) is less than the 
corotation radius, r CO i (the radius where the disc Keplerian fre¬ 
quency matches the neutron star spin frequency, which is 24 km 
for PSR J1023+0038). If the inner edge of the disc extends inside 
the pulsar’s light cylinder radius (ri c ; where the field lines attached 
to the neutron star synchronously rotate with the star at a velocity 
equal to the spee d of light) but lies outside r cor , ‘ propeller accre¬ 
tion' can occur ( llllarionov & Sunvaevlll975l ). In this case the pres¬ 
sure of infalling material is balanced by the magnetic field of the 
neutron star outside n c (r m »r co r) and the infalling matter does 
not reach the neutron star surface. The resulting centrifugal barrier 
created by the rotating neutron star expels the material from the 
system and the radio pulsar’s emission is qu enched l lSpruit & Taaml 
Il993l ; iRappaport. Fregeau. & Snruill2004l) . However, it is worth 
noting that the radiation pressure of a radio pulsar does not neces¬ 


sarily have to disrupt the disc and the disc can remain stable outside 
n c jEksi & Alparl20051) . 

The recent detection of X-ray pulsati ons in PSR J1023+0038 
during the accretio n active state l lArchibald et al.1 120141 ; 
[Bogdanov et al J 120141) implies that channeled accretion, simi¬ 
lar to that seen in higher-luminosity accreting millisecond X-ray 
pulsars, is occurring at a much lower accretion rate, implying 
that the inner edge of the accretion disc lies close to r co r and 
a propeller does not form. In this case material can accumulate 
near r COI and non-stationary accretion can occur as matter piles 
up around the intrinsically unstable magnetospheric boundary. 
Accretion discs accreting onto the magnetosphere of a rotating star 
can end up in a trapped state, in which the inner edge of the disc 
stays near r cor - The captured material can form a quasi-spherical 
shell jPringle & Reesl 1 19721: IShakura et all 1 20121) or a new disc 
struct ure, known as a ‘dead-disc’ a nd episodic accretion can 
occur ( ID’ Angelo & Spruitll201(]Ll2012l) . As noted by lPatruno et al] 
l i2014l) . variations in the mass-accretion rate of the accretion disc 
can move the inner disc radius by a factor of two or three. The 
viscous time-scale defines the time-scale for this drift and to 
reach a viscous time-scale of 10-100 s requires a region with an 
annulus of radius ~ 10-100 km. The thermal time-scale is given by 
(therm ~ (H/R) 2 ( v i s (where H and R are the height and radius 
of the disk) and for an inner thin disk with H/R <0.02 implies 
(therm <few seconds, much lower than the time-scales observed. 
Thus, in principle, fluctuations in the mass accretion rate can move 
r m outside r\ c ,which allows the radio pulsar to turns on, triggering 
a transition to the passive state, where the lower X-ray luminosity 
is produced by a shock between the pulsar wind and innermost 
accretion flow. An increase in the mass accretion rate pushes the 
inner edge of the disc back inside the light cylinder and turns off 
the radio pulsar. The X-ray pulsations are only observed in the 
active state and not in the passive state light curves, which suggests 
that the switching between the passive and active states results in 
transitions between a non-a ccreting pure propeller mode and an 
accreting trapped-disc mode jArchibald et alj2014i) . 

To explain the optical/X-ray timing and spectral properties of 
the passive and active state light curves, we propose a scenario in 
which the optical flux during the passive and active states arise from 
the viscously-heated and X-ray reprocessed disc, respectively, as 
well as from the irradiated secondary star. High energy photons 
(presumably X-rays) produced in the inner parts of the accretion 
disc, very close to the compact object, photoionize and heat the sur¬ 
rounding outer regions of gas in the accretion disc and secondary 
star. These regions later recombine and cool, producing lower en¬ 
ergy photons. This results in reprocessed optical and UV radiation 
arising from a volume of significant spatial extent, that is imprinted 
with the same variability as the X-ray signal. The steeper PDS at 
shorter wavelengths suggests more high frequency variability aris¬ 
ing from reagions closer to the compact object, as expected. An 
increase in the mass accretion rate pushes the inner edge of the disc 
inwards, past some critical radius which lies most likely near n c , 
and the radio pulsar turns off. During this accretion-powered active 
state, matter is accreted onto the neutron star surfac e via the mag- 
netic poles, resulting in X-ray pulsations ^Archibald et al.l2014l) . as 
well as irradiating the outer parts of the accretion disc. The buildup 
and release of mass in a trapped disc at some critical radius leads 
to clumpy accretion onto the neutron star’s surface, as depicted by 
the rectangular flat-bottom light curves. A decrease in the mass ac¬ 
cretion rate moves the inner edge of the accretion disc outwards 
allowing the radio pulsar to turn on, triggering a transition to the 
rotation-powered passive state, where the lower X-ray luminosity 
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is produced by a shock between the pulsar wind and innermost ac¬ 
cretion flow and the viscously-heated accretion disc. 


8 CONCLUSIONS 

Below we list the main results of this paper. 

(i) Our long-term time-resolved r'-band light curves taken be¬ 
tween 2014 January and April show an underlying sinusoidal mod¬ 
ulation due to the irradiated secondary star and accretion disc. We 
also observe superimposed rapid flaring on time-scales less than 
minutes with amplitudes of ~0.1-0.5 mag, and further additional 
large flare events on time-scales of ~5-60 min with amplitudes of 
~0.5-1.0 mag. The power density spectrum of the flare light curves 
is dominated by a red-noise component with a power-law index typ¬ 
ical of aperiodic activity in X-ray binaries. 

(ii) The UV and X-ray luminosities are strongly correlated sug¬ 
gesting that the same emission mechanism is powering part of the 
X-ray and UV emission. The power-law index of the X-ray/UV lu¬ 
minosity correlation is consistent with an irradiated accretion disc. 

(iii) For the large flare event observed on March 16 where we 
have colour information, we find that the event arises from a rel¬ 
atively optically thin region with temperature ~ 8700 K, Nn ~ 
6 x 10 22 nucleons cnk 2 and a radius of ~ 0.35 Rq, possibly from 
the accretion disc corona. 

(iv) On some nights we also observe sharp-edged, rectangular, 
flat-bottomed dips randomly distributed with orbital phase. These 
rectangular dips are similar to the mode-switching behaviour be¬ 
tween active and passive luminosity states observed in the X-ray 
light curves of other similar millisecond pulsars. The dips have a 
median duration and ingress/egress time of ~ 250 s and ~ 20 s 
respectively and we observe an increase in the blue colour ratio 
during the dips. This is the first time that the optical analogue of 
the X-ray mode-switching has been observed. 

To explain the timing and spectral properties of the active 
and passive state light curves, we propose a scenario in which 
the optical flux during the passive and active states arise from the 
viscously-heated and X-ray reprocessed disc and changes in the in¬ 
ner disc radius, which lies near r COI . An increase in the mass accre¬ 
tion rate leads to an accretion-powered active state, where matter is 
accreted onto the ne utron star surface via m agnetic poles, resulting 
in X-ray pulsations ( Archibald e t al.ll2014h . as well as irradiating 
the outer parts of the accretion disc. A decrease in the mass accre¬ 
tion rate moves the inner disc’s edge outwards, triggering a transi¬ 
tion to the rotation-powered passive state, where the lower X-ray 
luminosity is produced by a shock between the pulsar wind and the 
innermost accretion flow and the viscously-heated accretion disc. 
The optical/X-ray timing properties of the passive and active state 
light curves, can be explained in terms of transitions between a non¬ 
accreting pure propeller mode accretion and an accreting trapped- 
disc mode accretion with X-ray reprocessing, during which clumpy 
accretion from a trapped inner accretion disc near the corotation 
radius results in rectangular, flat-bottomed optical and X-ray light 
curves. 


ACKNOWLEDGEMENTS 

TS would like to thank Mallory Roberts and Rene Breton for useful 
discussions. This research has been supported by the Spanish Min¬ 
istry of Economy and Competitiveness (MINECO) under the grants 


AYA2010-18080, AYA2012-3 8700 and AYA2013-42627. PRG is 
supported by a Ramon y Cajal fellowship RYC-2010-05762. VSD, 
SPL and ULTRACAM are supported by STFC grant ST/J001589/1. 
Based on scheduled observations made with the William Herschel 
Telescope and the Isaac Newton Telescope operated on the island 
of La Palma by the Isaac Newton Group in the Spanish Observa- 
torio del Roque de Los Muchachos of the Instituto de Astrofisica 
de Canarias and on observations made with the IAC80 telescope 
operated on the island of Tenerife by the Instituto de Astrofisica de 
Canarias in the Spanish Observatorio del Teide. 


REFERENCES 

Alpar M. A., Cheng A. F., Ruderman M. A., Shaham J., 1982, Natur, 300, 
728 

Archibald A. M„ et al„ 2009, Sci, 324, 1411 

Archibald A. M., Kaspi V. M., Bogdanov S., Hessels J. W. T., Stairs I. H., 
Ransom S. M., McLaughlin M. A., 2010, ApJ, 722, 88 
Archibald A. M., Kaspi V. M., Hessels J. W. T., Stappers B., Janssen G., 
Lyne A., 2013, arXiv, arXiv:1311.5161 
Archibald A. M., et al, 2014, ApJ, arXiv, arXiv:1412.1306 
Amaud K. A., 1996, XSPEC: The First Ten Years, in Astronomical Data 
Analysis Software and Systems V, ed. Jacoby, G. H. and Barnes, J., 
ASPC, 101, 17 

Bogdanov S., Archibald A. M., Hessels J. W. T., Kaspi V. M„ Lorimer D., 
McLaughlin M. A., Ransom S. M., Stairs I. H., 2011, ApJ, 742, 97 
Bogdanov S., et al.. 2014, arXiv, arXiv:1412.5145 
Bond H. E„ White R. L.. Becker R. H„ O’Brien M. S., 2002, PASP, 114, 
1359 

Burrows D. N., et al., 2005, SSRv, 120, 165 

Campana S., Stella L., Kennea J. A., 2008, ApJ, 684, L99 

Cardelli J. A., Clayton G. C., Mathis J. S„ 1989, ApJ, 345, 245 

Church M. J.. Bafucinska-Church M„ 2004, MNRAS, 348, 955 

Coti Zelati F„ et al., 2014, MNRAS, 444, 1783 

D’Angelo C. R., Spruit H. C., 2010, MNRAS, 406, 1208 

D’Angelo C. R., Spruit H. C., 2012, MNRAS, 420, 416 

Deeg H. J., Doyle L. R., 2001, in Third Workshop on Photometry, ed. W. 

J. Borucki, & L. E. Lasher, 85 
Deller A. T„ et al., 2012, ApJ, 756, L25 
Dhillon V., et al., 2007, MNRAS, 378, 825 
Dhillon V. S., et al.. 2014, MNRAS, 444, 3504 
Ek§l K. Y.. Alpar M. A„ 2005, ApJ, 620, 390 
Ferrigno C„ et al., 2014, A&A, 567, A77 

Halpern J. P, Gaidos E., Sheffield A., Price-Whelan A. M., Bogdanov S., 
2013, ATel, 5514, 1 

Horne J. H„ Baliunas S. L„ 1986, ApJ, 302, 757 
Hynes R. I., et al., 2003, MNRAS, 345, 292 

Hynes R. I., Charles P. A., Casares J., Haswell C. A., Zurita C., Shahbaz 
T„ 2003, MNRAS, 340, 447 
Illarionov A. F., Sunyaev R. A., 1975, A&A, 39, 185 
Kerr G. S.. Fletcher L„ 2014, ApJ, 783, 98 
Kong A. K. H„ 2013, ATel, 5515, 1 
Lai D„ 2014, arXiv, arXiv: 1411.7995 

Linares M., Casares J., Rodriguez-Gil P., Shahbaz T., 2014, ATel, 5868, 1 

Linares M„ et al., 2014, MNRAS, 438, 251 

Linares M„ 2014, ApJ, 795, 72 

Linares M., et al., 2015, in prep. 

de Martino D., et al., 2013, A&A, 550, A89 

Papadakis I. E., Lawrence A., 1993, MNRAS, 261, 612 

Patruno A., Watts A. L„ 2012, arXiv, arXiv: 1206.2727 

Patruno A., et al., 2014, ApJ, 781, L3 

Pearson K. J., Horne K„ Skidmore W„ 2003, MNRAS, 338, 1067 

Predehl P„ Schmitt J. H. M. M„ 1995, A&A, 293, 889 

Press W. H., Teukolsky S. A., Vetterling W. T.. Flannery B. R, 1992 

Pringle J. E„ Rees M. J., 1972, A&A, 21, 1 

Rappaport S. A., Fregeau J. M., Spruit H., 2004, ApJ, 606, 436 


© 2014 RAS. MNRAS 000, ??-?? 




14 T. Shahbaz et al. 


Roming P. W. A., et al., 2005, SSRv, 120, 95 

Roberts M. S. E., 2011, Roberts, M. S. E. 2011, Radio Pulsars: An As- 
trophysical Key to Unlock the Secrets of the Universe (AIP Conf. Proc. 
1357), ed. M. Burgay et al. (Melville, NY: AIP), 127 
Russell D. M., Fender R. P., Hynes R. I., Brocksopp C.. Homan J., Jonker 
P. G., Buxton M. M„ 2006, MNRAS, 371, 1334 
Shahbaz T., Dhillon V. S., Marsh T. R., Zurita C.. Haswell C. A., Charles 
P. A., Hynes R. I., Casares J„ 2003, MNRAS, 346, 1116 
Shahbaz T., Russell D. M., Zurita C., Casares J., Corral-Santana J. M., 
Dhillon V. S., Marsh T. R„ 2013, MNRAS, 434, 2696 
Shahbaz T., Dhillon V. S., Marsh T. R., Casares J., Zurita C., Charles P. A., 
2010, MNRAS, 403, 2167 

Shahbaz T.. Dhillon V. S., Marsh T. R., Casares J., Zurita C., Charles P. A., 
Haswell C. A., Hynes R. I., 2005, MNRAS, 362, 975 
Shahbaz T., Hynes R. I., Charles P. A., Zurita C., Casares J., Haswell C. A., 
Araujo-Betancor S., Powell C., 2004, MNRAS, 354, 31 
Shakura N., Postnov K„ Kochetkova A., Hjalmarsdotter L., 2012, MN¬ 
RAS. 420, 216 

Spruit H. C., Taam R. E„ 1993, ApJ, 402, 593 
Stappers B. W„ et al., 2013, ATel, 5513, 1 
Szkody P, et al., 2003, AJ, 126, 1499 
Takata J., et al., 2014, ApJ, 785, 131 
Tauris T. M„ 2012, Sci, 335, 561 

Tendulkar S. P, et al., 2014, ApJ, arXiv, arXiv: 1406.7009 
Thorstensen J. R., Armstrong E., 2005, AJ, 130, 759 
Timmer J., Koenig M„ 1995, A&A, 300, 707 
Uzdensky D. A., 2004, Ap&SS, 292, 573 

Van Der Klis M„ 1989, H. Ogelman H., and E. P. van den Heuvel E.P. 
Timing Neutron Stars: proceedings of the NATO Advanced Study Insti¬ 
tute on Timing Neutron Stars, Kluwer Academic, New York, p.27 
van Paradijs J., McClintock J. E., 1994, A&A, 290, 133 
Vaughan S., 2005, A&A, 431, 391 
Walkowicz L. M„ et al., 2011, AJ, 141, 50 

White, N. E., Nagase, F., Parmar, A. N. , 1995, in X-Ray Binaries, W. H. 
G. Lewin, J. van Paradijs, E. P J. van den Heuvel (eds.) (Cambridge: 
Cambridge Univ. Press), 1 

Wilms J., Allen A.. McCray R„ 2000, ApJ, 542, 914 
Wijnands R., van der Klis M., 1998, Natur, 394, 344 
Zurita C., Casares J., Shahbaz T., 2003, ApJ. 582, 369 


© 2014 RAS, MNRAS 000, ??-?? 


